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Abstract 



We consider the effective potential in three-dimensional models with O(N) 
symmetry. For generic values of N, and in particular for the physically inter- 
esting cases N = 0, 1,2,3, we determine the six-point and eight-point renor- 
malized coupling constants which parametrize its small-field expansion. These 
estimates are obtained from the analysis of their e-expansion, taking into ac- 
count the exact results in one and zero dimensions, and, for the Ising model 
(i.e. N = 1), the accurate high-temperature estimates in two dimensions. 
They are compared with the available results from other approaches. We also 
obtain corresponding estimates for the two-dimensional O(iV) models. 
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I. INTRODUCTION 



The effective potential is widely used in the field-theoretic description of fundamental 
interactions and phase transitions. In field theory the effective potential is the generating 
functional of the one-particle irreducible correlation functions at zero external momenta. In 
statistical physics it represents the free-energy density T as a function of the order parameter, 
which, for spin models, is the magnetization M. Its global minimum determines the value of 
the order parameter, which characterizes the phase of the system. In the high-temperature or 
symmetric phase the minimum is unique with M = 0. As the temperature decreases below 
the critical value, the effective potential takes a double-well shape: the order parameter 
does not vanish anymore and the system is in the low-temperature or broken phaseQ The 
equation of state is closely related to the effective potential. It relates the magnetization M 
(i.e. the order parameter), the magnetic field H and the reduced temperature t oc T — T c . 
It is simply given by 



H 



dM' 



(1) 



In this paper we study the effective potential of O(iV) models. We recall that O(N) 
models describe many important critical phenomena in nature: liquid-vapour transitions 
in classical fluids, the helium superfluid transition, the critical properties of isotropic ferro- 
magnetic materials and long polymers. We will focus mainly on the small-renormalized-field 
expansion of the effective potential in the symmetric phase. Its coefficients are directly re- 
lated to the zero- momentum n-point renormalized coupling constants g n . The four-point 
coupling g = g A plays an important role in the field-theoretic perturbative expansion at 
fixed dimension ||, which provides accurate estimates of critical indices and universal ra- 
tios in the symmetric phase. In this approach any universal quantity is obtained from a 
series in powers of g ((^-expansion) which is then resummed and evaluated at the fixed-point 
value of g, g*. Accurate estimates of g* have been obtained by calculating the zero of the 
Callan-Symanzik /5-function associated with g (see e.g. Refs. PH91). These results have been 



substantially confirmed by computations using different approaches, such as e-expansion []T0 



high-temperature expansion (see, e.g., Refs. fIT|- |i4||K)|| ), Monte Carlo simulations (see, e.g., 
Refs. QjHni), l/N expansion 

Recently there has been a lot of interest in the problem of computing the higher-order cou- 
pling constants g 6 , g 8 , . . ., for the Ising model (see e.g. Refs. |^|5l00^ggH^M3i3)- 
Here we will study the issue for generic values of N. We will obtain satisfactory estimates of 
g 6 and g 8 , or, equivalently, of the ratios r 2 j = (fej/V -1 , from their e-expansion [^3|] (e = 4— d) 
within the 4 theory defined by the action 



S 



d d x 



]-d^(x)d^(j)(x) + ^m 2 <f) 2 + ^g ( 



<.2\2 



(2) 



Actually in the broken phase the double-well shape is not correct because the effective potential 
must always be convex. In this phase it should present a flat region around the origin. For a 
discussion see e.g. Refs. ]l]|§. 
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The e-expansion of the fixed-point value of r 2 j can be derived from the e-expansion of 
the equation of state, which is known to 0(e 3 ) 



i.e. 



N 



for the Ising model |2^J25[ 
and to 0(e 2 ) for generic values of N |26|, thus leading to series of the same length for 
r 2 j. Since the available series are short and have large coefficients increasing with j, their 
straightforward analysis does not provide reliable estimates, but gives only an indication of 
the order of magnitude. A considerable improvement can be achieved if one uses the results 
for d = 0, 1, 2, whenever they are available. This idea was employed in Ref. Effl to improve 



the estimates of the critical exponents of the Ising and self-avoiding walk models, and in 
Ref. |ZB| it was used to the study of the two-point function. In Ref. |TD[ it was generalized 
and successfully applied to the determination of the zero-momentum four-point renormalized 
coupling. In the present case the basic assumption is that the zero-momentum n-point 
renormalized couplings g 2 j, and therefore the ratios r 2 j, are analytic and quite smooth in 
the domain 4 > d > (thus < e < 4). This can be verified in the large- N limit. One may 
then perform a polynomial interpolation among the values of d where the constants r 2 j are 
known (d = 0, 1) or for which good estimates are available (d = 2 for N = 1, obtained from 
a high-temperature analysis), and then analyze the series of the difference. This procedure 
leads to more accurate estimates, which are consistent with those obtained by the direct 
analysis of the original e-series, but have a much smaller uncertainty. As a by-product of 
our analysis we also obtain relatively good estimates of r 6 and r 8 for two-dimensional 0(N) 
models. 

For N 1, most of the published results concern the renormalized four-point coupling 
constant g. As far as we know, estimates of g§ and g% have only been obtained from 
approximate solutions of the renormalization group equation [p|, and from the analysis 



. The latter results present a large uncertainty. The former 
which, in any case, are very difficult to estimate — and their 



of high-temperature series [[12 
are reported without errors — 
reliability is unclear. For instance, the estimates of g that are obtained using this method 
are in disagreement with the results of other computations. Therefore our new independent 
estimates of r§ and r§ for N ^ 1 represent our main results, and provide an important check 
of the above-mentioned calculations. 

For the Ising model we can compare our results with the estimates obtained using dif- 
ferent approaches: the (/-expansion at fixed dimension d = 3 pO|,P| that apparently provides 



the most precise results; a different analysis of the e-expansion based on the parametric 
representation of the equation of state ||; approximate solutions of the exact renormaliza- 
tion group equation [rgpi| ; high-temperature expansions [|TJ],[r^,|IB|] ; dimensional expansion 
around d = p^JT9|] ; Monte Carlo simulations [|1^,|17]]. Our final results are in good agree- 
ment with these estimates and their precision is comparable with that of the analysis of the 
(/-expansion. 

The paper is organized as follows. In Sec. |T| we introduce our notation and give some 
general formulae for the small-field expansion of the effective potential. In Sec. [TIT] we 
compute the effective potential in the large- N limit, the ratios r 2 j to 0(1/N), and give the 
exact values of r 6 and r§ for d — 1,0. Furthermore we present a high-temperature analysis 
of the two-dimensional Ising model, which provides accurate estimates of the first few r 2 j. In 
Sec. [TV] we present our analysis of the e-expansion of r 2 j. In Sec. |V] we compare our results 
with other approaches. In App. [A] we give some useful formulae relating the ratios r 2 j to 
the connected Green's functions. In App. |B| we report the e-expansion of r 2 j derived from 
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the e-expansion of the equation of state. In App. |C| some details of the calculations in one 
dimension are given. 



II. SMALL-FIELD EXPANSION OF THE EFFECTIVE POTENTIAL 

The free energy per unit volume can be expanded in powers of the renormalized magne- 
tization ip (i.e. the expectation value of the renormalized field <p r = 4>/\fZ): 

n<p) = ^ = ^(°) + l m V + ^r 4 ~w + >> 2 '- ;i J ' ! j^!ivr 21 . (3) 

where is the generating functional of one-particle irreducible correlation functions at 
zero external momenta, i.e. the effective potential of the renormalized theory. The mass 
scale m is the inverse of the second-moment correlation length, i.e. m = £ _1 and 

2 = 1 Jdx x 2 G{x) 
5 2d Jdx G(x) ' 1 ' 

where the function G(x) is defined by 

(0«(O)M*)> = ScfiG(x). (5) 

By rescaling (p as 

m (d-2)/2 

43 = («) 
in Eq. (|^), the free energy can be written as 

r(<p)-F(0) = —A(z) t (7) 
9 

where 
and 

In App. [A| we give some useful formulae to derive the constants r 2 j from the connected 
Green's functions. 

One can show that z oc t _/3 M, and that the equation of state can be written in the form 

H oc (10) 

oz 
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This relation can be exploited in order to derive A(z) from the equation of state, which is 
usually written in the form (see e.g. Ref. 0]) 

H = M 5 f(x) (11) 

where x = tM' 1 ^ . The function A(z) is thus given by 

^ = hoJf (x z-^) , (12) 

where the normalization constants h and Xq are fixed by the requirement that 

A{z) = \z 2 + ^ + 0{z% (13) 

The ratios r 2 j are obtained by expanding A(z) in powers of z. Notice that, since the function 
f(x) in Eq. ( ]TT| ) is regular at x = and nonzero, Eq. (|1^) implies A(z) ~ z s+1 for z — > oo. 

In the following we will be interested in the rescaled effective potential A(z) and we will 
calculate the fixed-point values of the first few coefficients r 2 j of its small- ,2 expansion. 

III. EXACT AND HIGH-TEMPERATURE RESULTS 

In order to get a qualitative idea of the properties of the effective potential, we consider 
its large-iV limit. It is easy to derive the large-iV limit of the rescaled effective potential 
A(z) from the corresponding equation of state p3f: 



H = M 5 (l + xf /{d ~ 2) , (14) 
where 6 = (d + 2)/(d- 2), x = tM~ 1 ^ and (3 = 1/2. One finds 

6 



^ d 



d-2 ^ d/(d ~ 2) 



l + —z 2 
. 12 



(15) 



from which 



*«- *t-AS M- 2L f)- (l6) 



i=l 



The constants r 2 j are (j — 2)th-order polynomials in e = 4 — d that have j — 2 real zeros at 
e = 2{i — with i = 1, — 2. Notice that in the limit j — > oo the zeros have e = 2 (i.e. 
d = 2) as an accumulation point. It is interesting to note the form of the large-iV limit of 
A(z) for integer d: 

A(z) = \z 2 + ^ 4 for d = 4, (17) 

A(z) = \z 2 + ±z 4 + ^z e for d = 3, (18) 

A(z) = 3 (V 2 / 6 - l) for d = 2, (19) 

A{z) = 6z 2 {12-z 2 f 1 for d=l. (20) 
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Notice that in the large- iV limit, for d = 3,4 (and in general for d = 2n/(n — 1) with 
integer n > 2) the effective potential is a polynomial in z, or equivalently in ip. In four 
dimensions only the first two terms are present so that the effective potential (in the variable 
z) coincides with the phenomenological expression of Ginzburg and Landau. Notice however 
that the rescaling (^) is not strictly defined in four dimensions since g — > in the critical 
limit. Therefore this simple expression is not valid in the original variable p and indeed 
logarithms of the magnetization appear in the four- dimensional effective potential |nj. In 



three dimensions also the ip term is present so that 5 = 5. As d — > 2 all terms become 
relevant. Of course this simple behaviour is peculiar of the large- N limit. For finite values 
of N all terms are present in the small-field expansion, and 5 can only be determined after 
resumming the series. 

One can also derive the 0(1/N) correction to Eq. ( TO ) from the corresponding 0(1/N) 
correction to the equation of state calculated^ in Ref. ||30|| . In d = 3 one obtains 



r6 = 6 



12.2556 „ / 1 



(21) 



67.3140 n / 1 \ 

r » = + \jp) • < 22 > 

1406.83 

rio = _ + (_), (23) 

etc.... For large values of j the coefficients of the 0(1/N) term in the 1/N expansion of T2j 
behave approximately as ~ (2j)!(— c) 3 ', where c is a constant: c ~ 0.40. The large coefficient 
of the 0(1/N) correction in r 6 indicates that the region where Eq. ( fUD may be a good 
approximation corresponds to very large values of N, say iV > 100. This is expected to be 
true also for r 2 j with larger values of j. 

The constants r 2 j can be computed exactly in one and zero dimensions. These results 
will be useful in our analysis of the e-expansion of r 2 j as explained in the introduction. 

Some details of the calculations for d = 1 can be found in App. y. Here we only give 
the results for tq and r§. For iV > 1 we have 

r _ 5 5N(N-m8N + 7) 

6 (iV + l)(iV + 4)(4iV- l) 2 ' 1 J 

175 _ 35N(N - l) 2 (256iV 3 + 3037A 7 " 2 + 1705A^ - 588) 
r§ ~ ~3 3(N + 1)(N + 4)(N + 6)(AN - l) 3 ' ( } 



For N < 1 instead 



r 6 = 5, (26) 
175 

rs = -z-. (27) 



2 We mention the presence of a misprint in the final expression of the 0(1/N) equation of state 
given in Eq. (29) of Ref. [ 30 1 : in the third term of the first line (27r) 2_<: / 2 should be replaced by 
2^ 2 ^/ 2 . 



6 



For N = 1 these expressions agree with the results of Ref. ||19|| . For N = oo they reproduce 
Eq. ([161) . For the Ising model we will also need the value of r 10 , which has been computed 
in Ref. 0: r 10 = 1225. 

For d = and N > 1, using the formulae reported in App. [A|, it is easy to obtainP] 

_ 10(iV + 8) 

r6 " 3(iV + 4)' (28) 

70(iV 2 + 14iV + 120) 

8 3(iV + 4)(7V + 6) ' 1 J 

280(10752 + 3136iV + 256iV 2 + 30iV 3 + iV 4 ) , , 

r io = TTt , ^9/ at- , g\ / Ar , o \ • ( 30 ) 



For N = 1 these results agree with those of Ref. ||19|| . It is not clear how to determine the 



value of T2j for N = 0. Unlike the case d = 1, we cannot prove that setting N = in the 
formulae obtained for iV > 1 provides the correct answer. 

For the two-dimensional Ising model reliable estimates of the first few r 2j - can be obtained 
from the analysis of their high-temperature expansion on the lattice. The basic reason 
is that the leading correction to scaling is analytic, since the subleading exponent A is 



expected to be larger than one |32| , |33f . Therefore the traditional methods of analysis of high- 
temperature series should work well. The series published in Refs. [|34"| , |35l for the lattice Ising 
model with nearest-neighbor interactions allow us to calculate r 2 j (more precisely, a high- 
temperature series whose value at the critical point is r 2 j) to 17th-order on the square lattice 
(for which (3 C = ArcTanh(\/2 — 1)) and to 14th-order on the triangular lattice (for which 



P c = ArcTanh(2 — v^))- In the analysis of these series we followed Ref. |H|, using several 
types of approximants, Pade, Dlog-Pade and first-order integral approximants. Table [I 
reports the results obtained on the square and on the triangular lattice. They are consistent 
with each other. Assuming universality, as final estimates for the two-dimensional Ising 
model we take 

r 6 = 3.678(2), (31) 
r 8 = 26.0(2), (32) 
rio = 275(15). (33) 

The error on the estimate of r 2 j increases with j and thus the analysis of the higher-order 
coefficients does not lead to reliable estimates. We mention that the high-temperature 
analysis of Ref. |l3j leaf] to r 6 = 3.679(8), which is perfectly consistent with our estimate. 



3 The calculation is easily done for the N- vector model: in this case one has a single field s with 
s ■ s = 1 and the Gibbs measure is simply ds 5(s- s — 1). 

4 Actually Ref. gives an estimate of Rq = x1/(X2Xg) = (10 — t^) 1 - 
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IV. ANALYSIS OF THE e-EXPANSION 



In this Section we will compute r 2 j for j = 3,4,5 using the e-expansion. The series in e 
of Tq, r 8 and r 10 are reported in App. [TJ. They were obtained from the e-expansion of the 
equation of state p6jp4|,p5| . Since the e-expansion is asymptotic, it requires a resummation 
to get estimates for d = 3, i.e. e = 1, which is usually performed assuming its Borel 
summability. 

The main point of our analysis is the use of the exact values of r 2 j for d = 0, 1 we have 
reported in the previous Section, and, for N = 1, of the precise two-dimensional estimates 
which have been obtained from the analysis of high-temperature series, see Eqs. ([TT]) , . 
(j33[). Indeed the constants r 2 j are expected to be analytic in the domain 4 > d > 0. This 
can be explicitly verified in the large-iV limit where the constants r 2 j are polynomials in 
e, cf. Eq. flip]). Moreover it was implicitly assumed in the dimensional expansion around 
d = done in Refs. [gg,]BJ. 

The idea of the method is the following: consider a generic observable and let R(e) be 
its expansion in e. Moreover suppose that the values of R are known for a set of dimensions 
ei,...,efc. In this case one may use as zeroth order approximation the value for e = 1 of the 
polynomial interpolation through e = 0, e±,...,ek and then use the series in e to compute the 
deviations. More precisely, let us suppose that exact values _R ex (ei), . . ., i? ex (efc) are known 
for the set of dimensions ei, . . ., e&, k > 2. Then define 

Q(e) = E 

and 

S{e)^-^—--Q{e), (35) 

and finally 

R [mp (e) = [Q(e) + S(e)]]l(e-e l ). (36) 

8=1 

One can easily verify that the expression 

[Q(e)+S(0)]n(e-e,) (37) 

i=l 

is the fc-order polynomial interpolation through the points e = 0, e\, e&. The resummation 
procedure is applied to S(e) and the final estimate is obtained by computing i?i mp (e = 1). 

Since the series of r 2 j begins with a term of order e, we analyze the quantity r 2 j/e. 
Notice that, as a consequence, the interpolation formula fl3"7|) actually uses the value of the 
derivative of r 2 j in four dimensions. If the interpolation is a good approximation one should 
find that the series which gives the deviations has smaller coefficients than the original one. 
Consequently one expects that also the errors in the resummation are reduced. We find 



R P 



n 



e,: - U 



-II j=l,j& 



(34) 
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that, as expected, the coefficients of the corresponding series S(e) decrease in size with k, 
the number of exact values that are used to constrain the series. This fact was also shown 
in Ref. |T(| for the case of the four-point renormalized coupling. 

provide further support to our constrained analysis. 



The large- N results of Eqs. 
Indeed one may consider the simple polynomial interpolation (which uses the values of r^j 



in d = 0, 1 and the value of its derivative in d — 4) evaluated at d 
large- N expansion with the exact one. One finds (for N > 1) 



3, rJy , and compare its 



5(16iV 5 + 402iV 4 + 1734iV 3 + 539iV 2 - 669iV + 84) 



5 
6 



Q(N + 1)(N + 4)(iV + 8)(4JV - If 
+ 8N N 2+ \N 3 , 



and 



„int 



75.65 988.4 / 1 



(38) 



(39) 



Comparing Eq. (^) with Eq. (pTD, one sees that r™* gives the exact result for N = oo. 
Moreover, also the 0(1/N) correction is closely reproduced: indeed the coefficient of the 
0(1/N) term in rg nt is 10.52 to be compared with the exact value 10.21 of Eq. (pip. Therefore 
in the large- N limit 7g nt provides an estimate of r$ with a relative error which behaves as 



Eq. 
for 



0.37/iV: Tg nt becomes increasingly accurate as iV — > oo. The same discussion applies to 
the coefficient of the 1/N term in Eq. (^) is very close to the exact one —67.31, cf. 
(f22|). Therefore also in this case the interpolation rg nt provides good estimates of r 8 : 
N — >• oo the relative error is 12%. 

The analysis of the series S(e), cf. Eq. (|35D , can be performed by using the method 
proposed in Ref. M, which is based on the knowledge of the large-order behaviour of the 
series. It is indeed known that the n-th coefficient of the series behave as ~ (— a) n r(n+6 + l) 
for large n. The constant a, which characterizes the singularity of the Borel transform does 
not depend on the specific observable; it is given by |36| , |37| a = 3/(N + 8). The coefficient 
bo depends instead on the series one considers. Given a quantity R with series 



R(e) 



fc=0 



(40) 



we have generated new series R p (a, b; e) according to 



f poo 

RJa, b;e) = Y Bda, b) / dt t b e~ l 

t^n JO 



u(et) k 



where 



fc=0 



u\x) 



[1 - u(et)] a 



\/l + ax — 1 
\/l + ax + 1 



(41) 



(42) 



The coefficients Bk(a,b) are determined by the requirement that the expansion in e of 
R p (a, b; e) coincides with the original series. For each a, b and p an estimate of R is simply 
given by R p (a, b;e = l). 
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For the Ising model, where the available series are of order 0(e 3 ), we followed Ref. 
in order to derive the estimates and their uncertainty. We determine an integer value of b, 
b opt , such that 



Rs{oi, b opt ; e = 1) « R 2 (a, b opt ; e = 1) (43) 

for a < 1. b opt is the value of b such that the estimate from the series to order 0(e 3 ) is 
essentially identical to the estimate from the series to order 0(e 2 ). In a somewhat arbitrary 
way we have then considered as our final estimate the average of R p (a, b;e = 1) with 
— 1 < a < 1 and —2 + b opt < b < 2 + b opt . The error we report is the variance of the 
values of R^a.b^e = 1) with — 1 < a < 1 and [b opt /3 — lj < b < |~4& opt /3 + 1]. This 
procedure is ad hoc, but provides estimates that are all consistent among each other. In 
order to test the method, in Ref. |T(| , the procedure was applied to the determination of the 
critical indices and it provided estimates and error bars in substantial agreement with the 
results of other authors. Therefore we believe that our error bars are reasonable, although 
one should be cautious in giving them the standard statistical meaning. 

The results of our analysis for the Ising model, corresponding to N = 1, are presented 
in Table y. We report various estimates of r 6 , r 8 and r 10 obtained from an unconstrained 
analysis and constrained analyses in various dimensions. They are all consistent. As ex- 
pected, the error decreases when additional lower dimensional values are used to constrain 
the analysis: the error of the unconstrained analysis is approximately an order of magnitude 
larger than the error of our best result that uses the known values at d — 0, 1, 2. In Figs. [I] 
and |2] we show respectively r 6 and r§ as a function of d. There we plot the polynomial 
interpolations through the known values of r^/e at d — 4, 2, 1, 0, and the results of our con- 
strained (d = 2, 1, 0) e-expansion analysis. Their comparison shows how well the polynomial 
interpolation works. 

We have also repeated the analysis in two dimensions. In this case, of course, it is 
more difficult to get precise estimates: the unconstrained expansion gives results with large 
errors and it is therefore practically useless. Better estimates are obtained constraining 
the expansion in one and zero dimensions. The results for these two cases are reported in 
Table |I[ The final estimates for r$ and r§ are in very good agreement with the much more 
precise results obtained from the high-temperature analysis. The result for rio is instead 
significantly lower than the high-temperature estimate. It should be noted however that the 
series in e for r 10 has very large coefficients and the estimates show large fluctuations with 
the parameters b and a. Therefore it is not clear if our algorithm to determine the error 
bars is working properly here. For this reason, in this case we believe the high-temperature 
estimate to be more reliable than the e-expansion result. 

For generic values of N, the series is one order shorter and we have only two non-trivial 
terms. The procedure we presented above cannot be applied and we used a different method. 
For each value of b, averaging over — 1 < a < 1, we obtain three estimates: the first one 
is the result of the analysis of the unconstrained series, the second and the third one the 
results from the series constrained respectively in d — 1 and in d — 1, 0. We then compute 
the weighted average and the corresponding x 2 . The optimal value b opt is chosen as the 
value of b with the smallest x 2 ■ Once b opt has been chosen, we calculated mean values and 
relative spreads varying b and a according to the algorithm we presented above. Again, we 
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must stress that the procedure is completely arbitrary, and, since the series has only two 
terms, one should be very cautious in interpreting the spread of the approximants as an 
error bar. Our three-dimensional results for r 6 and r 8 are reported in Table [TIT]. There we 
report the values of the simple polynomial interpolation r^f (see Eq. (j38[) for rg nt ), and the 



results from an unconstrained analysis and constrained analyses in various dimensions (the 
errors reported are the spreads of the approximants varying b and a). The corresponding 
results for d = 2 are presented in Table |TV| . 

To understand the reliability of our procedure we applied the above analysis to the case 
N — 1, i.e. to its 0(e 2 ) series of r 2 j constrained at d — 0, 1, and compared the results^ (see 
Table PTI| ) with the previous ones reported in Table [II]. In three dimensions the estimates 
of r& and r 8 that have been obtained using this method are respectively lower (higher) by 
half an error bar (resp. one error bar) than the corresponding estimates obtained from the 
longer series. So we expect the three-dimensional results for r 8 to be somewhat higher than 
the true result, at least for small values of N. In two dimensions there is good agreement 
for r 6 while r 8 differs by one error bar. Therefore, for AT = 1, it seems that our procedure 
gives reasonable estimates and error bars. 

A rigorous check of the error bars can be done in the large- A" limit. In the case of r^, 
consider the analysis which provides the best estimates, the constrained one in d = 1,0. 
Here r 6 is estimated from the expansion 



re _ 5 
e ~ 6 



1 10.63 + 0.37e 2 1 /Ar2 

1 + + 0(e 2 , l/N 2 



(44) 



This equation shows that inclusion of the 0(e) term gives an estimate of which is larger 
than Tg 11 *, in contrast with what we expect on the basis of the exact result (pT|). Clearly the 
deviation r 6 — rg nt is positive for small values of e and negative for e = 1. Such a behaviour 
cannot be reproduced by a single term in e and thus, at least for large values of N, the 
inclusion of the 0(e) in Eq. fl4~4|) worsens the final estimate. As a consequence the spread of 
the approximants largely underestimates the correct error. The previous discussion shows 
that this is not the case for N = 1 and thus we expect the problem to appear for some 
intermediate value of A" that is unfortunately unknown. We have decided to be conservative 
and we have simply assumed that for all values of A" the spread of the approximants is an 
unreliable estimate of the error. Thus, in Table [V| that summarizes our final results, the error 
we report for r 6 is the difference between the final result and the value of the interpolation. 
This method provides a reasonable estimate of the uncertainty in the large- A" limit, but, 
for the reasons we mentioned above, we believe it overestimates the errors in the small- A" 
region (this is evident in the N = 1 case). Of course, the two-dimensional estimates of r 6 
may have the same problem. In this case however the spread of the approximants is larger 
and includes r 6 — r 6 nt . Therefore, whatever method we use, we get essentially the same error 
bars. 



5 For iV = 1 we have also applied this method using the two-dimensional results. In this case 
we consider six estimates, from the analyses of the series quoted in Table |IJ From the series to 
order 0(e 2 ) (resp. 0(e 3 )) we get r 6 = 2.048(22) (resp. r 6 = 2.061(14)) and r 8 = 2.9(7) (resp. 
r 8 = 2.3(5)). 
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The same analysis can be repeated for r§, which is estimated from 

-75.65 + 6.35. +0(t2 _ 1/jy2) (45) 

In this case the 0(e) term has the correct sign. For e = 1 it provides estimates that differ 
from the correct result, Eq. (|22]), by 2%, which is indeed the size of the error bar of our final 
estimates. Therefore the error bars obtained for r§ are correct for large values of N. Since 
they are also reasonable for iV = 1 we believe them to be reliable for all values of N. 



V. CONCLUSIONS AND COMPARISON WITH OTHER CALCULATIONS 

We have studied the small-renormalized-field expansion of the fixed-point effective po- 
tential in the symmetric phase of three-dimensional O(N) models. The coefficients of this 
expansion are related to the zero- momentum n-point renormalized couplings g n . By prop- 
erly rescaling the effective potential one may re-express the small-field expansion in terms 
of the ratios r 2 j = gij/g^ 1 with j > 3, cf. Eq. ([8]). We have derived the e-expansion of r 2 j 
from the e-expansion of the equation of state which has been computed to order 0(e 3 ) for 
the Ising model [f?|, and to 0(e 2 ) for iV ^ 1 When, as in the case of the ratios r 2 j, the 



quantity at hand is expected to be analytic and smooth as a function of d, one may use exact 
results (or precise estimates) for lower-dimensional O(N) models in order to constrain the 
analysis of the corresponding e-expansion. For this reason we have computed the first few 
r 2 j exactly in one and zero dimensions for all values of N and we have estimated the same 
quantities in two dimensions for N = 1 from the high-temperature series of Refs. |34|,|35 



The constrained analyses of the available e-series of r 2 j, according to the procedure outlined 
in Sec. [TV], allowed us to achieve a considerable improvement with respect to their standard 
resummation, and led to satisfactory estimates of the first few r 2 j. Using the accurate esti- 
mates of the fixed-point value of g = g% which can be found in the literature (see Sec. |), one 
can extract the fixed-point value of the zero-momentum 2j-point renormalized couplings g 2 j 
from the relation g 2 j = g-*~ l r 2 j. 

Let us compare our results with the available estimates from other approaches. For 
N ^ 1 there are not many published results: we are only aware of the estimates of g±, g$ 
and gs presented in Refs. |jl8| , [T2|j . Table |V| presents a summary of all the available (as far as 
we know) estimates of and rs for several values of N ^ 1. Ref. |l8j uses a renormalizat ion- 



group approach in which the exact RG equation is approximately solved (no estimates of the 



errors are presented there). Ref. |12| instead derives g^ and g§ from their high-temperature 
expansion in the lattice iV-vector model. These estimates are in reasonable agreement with 
our results. One can also compare the results of our analysis for large values of iV with the 
1/N expansion to 0(1/N). This comparison shows a substantial consistency, although the 
region where the 0(1/ N) approximation of r 2 j is effective corresponds to very large values 
of N, say N > 100. 

We have also computed r 6 and r 8 in two dimensions. We recall that the two-dimensional 
N- vector model is asymptotically free for N > 3, and has a Kosterlitz-Thouless transition 
for N = 2. 

For N = 1, the Ising model, many works have been devoted to the study of the effective 
potential, exploiting various approaches. Table |VT] presents a summary of all the available 
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(as far as we know) estimates^ of r 6 , r 8 and r w . Apparently the most precise results are 
those obtained in Ref. ||. They have been derived from the analysis of the g-expansion 
at fixed dimension of the effective potential calculated to five loops p8| , |39| (in Table [VT] 
we refer to this approach by d = 3 g-exp.). Ref. || presents also an analysis of the 0(e 3 ) 
expansion using the parameteric represention of the equation of state (in Table [VT] we refer 
to this approach by e-exp. PREQ). In this case the e-expansion is used to estimate the coef- 
ficients of the expansion of the function^ h(9) characterizing the parametric representation 
of the equation of state (see e.g. Ref. [|]). Then the constants F 2 j-i = r 2j /(2j — 1)! are 
obtained using their relation to the expansion of h(8). Precise estimates of r 2 j have also 
been obtained in Ref. |HJ (see also Ref. by approximately solving the exact renormal- 
ization group equation (ERG), although the estimate of g by the same method is not equally 
good. Additional results have been obtained from high-temperature expansions |T^,|T2|jn 



and Monte Carlo simulations of the lattice Ising model. The high-temperature re- 

sults are in substantial agreement with the field-theoretic estimates. The apparent small 
discrepancy of the estimates of r 6 of Refs. [|T3| , p2| , which were obtained by using different 
lattice formulations of the Ising model, is probably due to the presence of confluent singu- 
larities which are not properly handled by standard approximants. The results of Ref. |14j 



come from an analysis of 0(f3 17 ) series on the cubic lattice, and they have been obtained by 
using the Roskies transform ^(| and suitably biased integral approximants [5T] which take 
into account the leading confluent singularity. High-temperature techniques have also been 
used to obtain a dimensional expansion around d = (d-exp.) of the Green's functions. 
The analysis of these series provides estimates of #4 and g§ ||19|| . The corresponding value 
of r§ is however smaller than the field-theoretic estimates, although the value found for g± 
is in substantial agreement []. The Monte Carlo results do not agree with the results of 
other approaches, especially those of Ref. JT7]]. But one should consider the difficulty of such 
calculations on the lattice due to the subtractions that must be performed to compute the 
irreducible correlation functions. In Ref. \I5\ estimates of #4 and g$ are obtained by looking 
at the probability distribution of the average magnetization. The discrepancy in this case 
may come from the 0(</> 6 ) polynomial approximation of the potential used to fit the data, 
or from possible finite-size effects. 

Our results for the Ising model are much more precise than those obtained for N ^ 1. 
This is due essentially due to two reasons: one additional order is known in the e-series 



When the original reference reports only estimates of g 2 j (see Refs. [12,19, 15| JT7|1 ) , the errors we 
quote for r 2 j have been calculated by considering the estimates of g 2 j as uncorrelated. 

7 The coefficients of the expansion of h(9) were estimated by setting e = 1 in their e-series. The 
reported errors take into account the uncertainty on the critical exponents 7 and 0, and the error 
on the normalization parameter p that relates r 2 j and the low-9 expansion of h(9). Other sources 
of error are neglected and therefore the final error may be underestimated. 

8 From the values of 54 and g% reported for the two dimensional model one derives rg = 3.12(12), 
to be compared with our strong-coupling result r§ = 3.678(2). 
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of r 2 j] beside exact results for d = 0,1, good estimates have been obtained in d = 2 by 
an analysis of the available high-temperature series, which can be used in the constrained 
analysis of the e-expansion. The precision of our results is comparable with that obtained 
from the analysis of the ^-expansion. Moreover there is a substantial consistency among the 
various approaches whose results are reported in Table |VJ. 



APPENDIX A: 



In this appendix we give some useful formulas relating r 2 j to the connected Green's 
functions evaluated at zero momentum. The coefficients of the expansion of A(z) around 
z = can be written in terms of one-particle irreducible correlation functions 

r 2 ^ = r£'2,.. a ,. a ,.(o,...,o), (Ai) 



as 



s = (W (iv + 2)^ r^rT (A2) 



gj-i 2^- 1 j!n£ 2 1 (iV + 2i) r 4 
In terms of the zero-momentum connected Green's functions 



X2j = ( s a 1 (0)s ai (x 2 )...s a] (x 2j ^ 1 )s a] (x 2j )) c , (A3) 

X2,...,X2j 



one then has 



Q _5(iV + 2)^ 

3(A + 4) xl V ' 

oon 280(A + 2) X6 X2 , 35(A + 2) 2 X sXl rA ,s 
r8 " 28U 3(A + 4) xl + 9(A + 4)(A + 6) xl ' { } 

7700(A + 2) X eX2 , 350( N + 2) 2 X \x\ 

'10 — 104UU — — r k h 



(A + 4) xl (A + 4) 2 xl 
1400(A + 2) 2 xsxl 35(A + 2) 3 X ioxl 



3(A + 4)(A + 6) xl 3(A + 4)(A + 6)(A + 8) x\ 



(A6) 



APPENDIX B: 



Here we present the e-series of r 2 j we used in our analysis. The e-expansion of r 2 j can be 
obtained from the equation of state which is known to 0(e 3 ) for the Ising model [24|,25,9|, 
and to 0(e 2 ) for generic values of N |26| . 



For generic values of N one finds 

5(A + 26) 4.37395A 2 + 55.6177A + 615.008 2 _ „ 



6(A + 8) (A + 

35(A + 80) 35(A 3 + 67.6582A 2 + 1661.61A+ 11634.7) 2 ^, SN 

-e+ — — — — — -e A + 0(e 6 ), (B2) 



18(A + 8) 18(A 
35(A + 242) 245(A 3 + 108.389A 2 + 4780.44A + 35830.0) 2 , N 
no = 3(A + 8) £ 12(A + 8)3 6 + 0(6 } - (B3) 



14 



For the Ising model, where one additional order is known, we have 



where 



5 25 2 /20C(3) 5A 310^ 



27 



V 9 



9 729 



0(e 4 



35 1925^ _ p50C(3) _ 35A + 18655^ ^ + 4 



— e + 
2 54 

13685 



V 9 



1458 
406945 \ 



r 10 = 315e - —^-e 2 + (^1260C(3) - 420A + ~^-) e 6 + 0(e 



r V3 



2tt 2 



1.17195 



and C(3) « 1.20206. 



(B4) 
(B5) 
(B6) 

(B7) 



APPENDIX C: 

In this appendix we present some details of the computation of r 6 and r 8 in d = 1. We 



use the formalism of Ref. [|42] which is based on the expansion of the Boltzmann weight in 
hyperspherical harmonics. Explicitly, given a generic nearest-neighbour Hamiltonian 



where s r ■ s n 



H = - ^ ^(4 • s x+ i) 

X 

1, we expand the Boltzmann weight as 



1 + H v N , k (P)Y N>k (si) ■ Y N , k (s 2 ) 



k=l 



(CI) 



(C2) 



where Y/v,fc(s) are the 0(A) hyperspherical harmonics. The coefficients vn,u{P) depend on 
(3 and on the explicit form of the interaction V(x). For the standard nearest-neighbour 
interaction, V(x) = x and 



VN,k(P) 



l N/2+k-l 



lN/2-x(P) 

where I n (/3) is a modified Bessel function. 

A rather lengthy computation gives the following results: 

12(A- 1) 



(C3) 



5-4 = 6- 



(A + 2) 



-R2 



g 6 = 180 + 180 



g 8 = 12600 + 2520 



(N-l) 
(A + 2) 2 
(A- 



-5(N + 2)R 2 



3A 2 
A + 4 



(A + 2) 3 



R\R ?J + (7A - A)R\ 
-44(A + 2) 2 i? 2 - 54 iV2 A ( r iV + 2) ^3 



(C4) 
(C5) 



-24 



+12 



N 3 {N+ l)(A + 2) 
(A + 4) 2 (A + 6) 

A + 4 2 



A + 4 

RjRjR 4 + Q N2 [ N + A 2) R\R\ + 8(A + 2)(17A - 11)^ 



A 4 



(A + 4) : 



A + 4 

-RlRl - 2(71A 2 - 82A + 20)i^ 



(C6) 
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Here Rj is the critical value of the ratio mi/rrij where rrij is the mass in the spin-j channel 
defined by the large-x behaviour of the spin-j correlation function, 



~ rrij \x\ 



(Ynj(so) ■ Y NJ (s x )} 
for \x\ — > oo. For the standard nearest-neighbour interaction we have 

N — 1 



(C7) 



R i 



j(N + j-2) 




for > 1 
for N < 1 



(C8) 



Substituting in the previous formulae, we get the results reported in Sec. [Til]. 
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TABLES 



TABLE I. Two-dimensional Ising model. We report the estimates of re, r$ and rio obtained 
from the analysis of the 17th order strong-coupling series on the square lattice and 14th order series 
on the triangular lattice (HT), and from the analyses of the e-expansion constrained at d = 1 and 



at d = 0, 1 (e-exp.). 





HT 




e-exp. 






square 


triangular 


d = 1 


d = 0,1 


re 


3.677(2) 


3.678(1) 


3.67(9) 


3.69(4) 


rs 


25.9(5) 


26.0(1) 


24.2(2.2) 


26.4(1.0) 


Ho 


269(11) 


279(11) 


131(81) 


171(31) 



TABLE II. Three-dimensional Ising model. Estimates of re,r§, and r\Q from the polynomial 
interpolation ("hit"), from an unconstrained analysis of the e-expansion, "unc", and constrained 
analyses in various dimensions. For the analyses which use the estimates in d = 2 we report two 
errors: the first one gives the uncertainty of the resummation of the series, the second one expresses 



the change in the estimate when the two-dimensional result varies within one error bar. 





hit 


unc 


d = 1 


d = 0, 1 


d = 2 


d= 1,2 


d = 0,1,2 


re 


2.092 


2.106(79) 


2.058(35) 


2.063(24) 


2.059(25+0) 


2.060(15+1) 


2.058(11+1) 




1.31 


0.4(2.4) 


1.93(89) 


2.65(63) 


2.23(60+2) 


2.53(39+3) 


2.48(27+5) 


'Id 


35 


-98(120) 


-7(67) 


15(38) 


-7(41+2) 


-8(18+5) 


-20(13+7) 



TABLE III. Three-dimensional estimates of r§ and r$ for various values of N from the polyno- 
mial interpolation ("int"), cf. Eq. (^) for tq, from an unconstrained analysis of the e-expansion 
and constrained analyses in various dimensions. 



N re rg 





int 


unc 


d = 1 


d = 0, 1 


int 


unc 


d = 1 


d = 0, 1 





2.42 


2.08(22) 


2.07(12) 




-7.8 


8(11) 


6.4(5.0) 




1 


2.17 


2.03(17) 


2.01(9) 


2.03(6) 


-0.6 


5.2(8.7) 


4.8 (3.6) 


4.7(1.9) 


2 


2.05 


1.94(14) 


1.91(6) 


1.94(5) 


-0.4 


3.6(7.0) 


3.2 (2.8) 


3.5(1.3) 


3 


1.93 


1.86(11) 


1.82(5) 


1.84(4) 


-1.2 


2.4(5.8) 


2.0 (2.2) 


2.1(1.0) 


4 


1.82 


1.76(9) 


1.74(4) 


1.75(3) 


-1.7 


1.5(4.9) 


1.1 (1.9) 


1.2(1.0) 


8 


1.500 


1.525(55) 


1.516(17) 


1.517(8) 


-2.5 


-0.7(2.7) 


-0.72(94) 


-0.74(44) 


16 


1.296 


1.304(26) 


1.296(6) 


1.291(1) 


-2.3 


-1.3(1.4) 


-1.29(39) 


-1.37(18) 


32 


1.104 


1.115(11) 


1.112(2) 


1.108(1) 


-1.6 


-1.4(7) 


-1.18(13) 


-1.22(6) 


18 


1.025 


1.036(6) 


1.032(1) 


1.029(1) 


-1.2 


-1.2(5) 


-0.96(7) 


-0.98(3) 
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TABLE IV. Two-dimensional estimates of rg and for various values of N from the polynomial 
interpolation of the known values at d = 4, 1, 0, from analyses constrained at d = 1 and at d = 0, 1. 



In the large- N limit r§ = 5/3 and r$ = 35/9. 



N 




re 












hit 


d=l 


d = 0,1 


int 


d=l 


d = 0, 1 





3.87 


3.7(9) 




10.5 


33(11) 




1 


3.78 


3.69(26) 


3.70(10) 


23.5 


28.7(7.5) 


28.6(2.6) 


2 


3.60 


3.51(22) 


3.54(7) 


20.9 


24.9(6.5) 


25.1(2.0) 


3 


3.38 


3.32(19) 


3.33(6) 


16.7 


20.3(5.4) 


20.3(1.7) 


1 


3.20 


3.15(15) 


3.15(5) 


13.6 


17.0(4.5) 


16.8(1.4) 


8 


2.73 


2.70(9) 


2.71(3) 


7.5 


9.3(2.2) 


9.3(7) 


16 


2.330 


2.34(4) 


2.325(5) 


4.3 


5.1(1.0) 


5.2(3) 


32 


2.045 


2.06(2) 


2.049(3) 


3.4 


3.82(40) 


3.81(12) 


48 


1.932 


1.95(1) 


1.936(2) 


3.3 


3.65(22) 


3.60(7) 



TABLE V. Summary of the available estimates of rg and rs for several values of N in three 
dimensions. 



N 


e-exp. 


re 

ERG |H HT @ 


1/N 


e-exp. 


^8 

ERG II 


1/JV 





2.1(3) 






6(5) 






2 


1.94(11) 


1.83 2.2(6) 




3.5(1.3) 


1.45 




3 


1.84(9) 


1.74 2.1(6) 




2.1(1.0) 


0.84 




4 


1.75(7) 


1.65 1.9(6) 




1.2(1.0) 


0.33 




8 


1.52(2) 




2.11 


-0.7(5) 




-8.41 


16 


1.291(5) 




1.47 


-1.4(2) 




-4.21 


32 


1.108(4) 




1.15 


-1.22(6) 




-2.10 


48 


1.029(4) 




1.046 


-0.98(3) 




-1.40 


100 


0.934(3) 


0.89 


0.9355 


-0.575(9) 


-0.64 


-0.67 
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TABLE VI. Summary of the available estimates of re, r% and rio for the three-dimensional 
Ising model. The error we quote for the result of our e-expansion analysis has been calculated by 
considering the two errors reported in the last row of Table |fj as uncorrelated. 



e-exp. [this paper] 
e-exp. PREQ § 
d = 3 5-exp. H 



ERG [21 



ERG |8| 



HT p| 



HT |3| 



HT 22 



HT g2| 

d-exp 



m 



MC p| 



MC |7| 



2.058(11) 

2.11(5) 

2.054(7) 

2.064(36) 

1.94 

1.99(6) 

2.157(18) 

2.25(9) 

2.5(5) 

1.54(26) 

2.72(23) 

3.26(25) 



?'8 



2.48(28) 

2.27(15) 

2.50(25) 

2.47(5) 

2.18 

2.7(4) 



12(2) 



-20(15) 
-11.6(7) 
-22(15) 
-18(4) 

-4(2) 
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FIGURES 



FIG. 1. For the Ising model we plot r§ as a function of d. The continuous line represents the 
polynomial interpolation through the known results at d = 4, 2, 1, 0. The bars are the results (with 
their uncertainty) of our constrained (d = 2, 1, 0) e-expansion analysis. 

FIG. 2. For the Ising model we plot r$ as a function of d. The continuous line represents the 
polynomial interpolation through the known results at d = 4, 2, 1, 0. The bars are the results (with 
their uncertainty) of our constrained (d = 2, 1, 0) e-expansion analysis. 
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